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Abstract 

Maximum likelihood methods, while widely used, may be non-robust due to disagreement between the 
assumptions upon which the models are based and the true density probability distribution of observed data. 
Because the Empirical Characteristic Function (ECF) is the Fourier transform of the empirical distribution 
function, it retains all the information in the sample but can overcome difficulties arising from the likelihood. 
This paper discusses, the ECF method proposed by Knight and Yu (2002), to estimate simultaneously the 
parameters for stable seasonal fractional ARIMA processes. Under some assumptions, we show that the 
resulting estimators are consistent and asymptotically normally distributed. For comparison purpose, we 
consider also a Two-Step Method (TSM) including in the first step the MCMC Whittle method developed 
by Ndongo et al. (2010), and in the second the MLE method introduced by Alvarez and Olivares (2005), 
to estimate the innovation parameters. The performance of the two methods is discussed under different 
parameter settings, using a Monte Carlo simulation. 

Keywords: Seasonal Fractional ARIMA, Stable distributions, ECF estimate, Whittle estimate, Markov 
Chains Monte Carlo, Two-Step Method. 



1. Introduction 

Seasonal Fractional ARIMA time series with stable innovations were introduced by Diongue et al. UK 
Since it is a direct generalization of the fractional ARIMA model of Kokoszka and Taqqu 1 15], it allows 
one to take into account three stylized facts: long range dependence, seasonality and high variability, often 
encountered in financial data. The use of processes with infinite variance has received a great deal of inter- 
est in the statistical literature as noticed in Brockwell and Davis [31]. Examples where such models appear 
to be appropriate have been found by Stuck and Kleiner 02511 . who considered telephone signals, and Fama 
, who modeled stock market prices. 

To solve the parameter estimation problem of the stable ARFISMA processes, Ndongo et al. lfl6ll propose 



several estimation procedures: the semiparametric method developed by Reisen et al. 12 1H . the classical 
Whittle estimate and the MCMC Whittle method which is based on the evaluation of the Whittle likelihood 
function using the Markov Chains Monte Carlo method (MCMC, in short). They study the behavior of these 
methods through Monte Carlo simulations, and the results show, in general, that the MCMC Whittle method 
is the best one. Unfortunately, this method only estimates the long-memory and short memory parameters. 
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Moreover, for inference purposes, estimation of the innovation parameter is required. The usual response to 
such difficulties is that the density functions of the processes cannot be written in a closed form in the sense 
that it is not expressible in terms of known elementary functions. As a result, ML estimation is often very 
difficult (for detailed discussion, see Calder and Davis [4]). Estimation of the stable ARFISMA model via 
alternative methods such as Quasi-Maximum Likelihood (QML) method and Generalized Method of Mo- 
ments (GMM) also presents difficulties. For example, QML is infeasible because the variance of the error 
term may be infinite. For GMM care must be taken when choosing moment conditions because the stable 
distribution does not have a finite absolute moment of order higher than the tail index a. Consequently, 
another alternative is to use the Empirical Characteristic Function (ECF) method. 

A main objective of this paper is to estimate simultaneously the parameters of the symmetric cr-stable 
ARFISMA processes using the Empirical Characteristic Function (ECF) method. The asymptotic proper- 
ties of the ECF estimators are established under some assumptions. Monte Carlo simulations are performed 
to study the finite sample properties of the ECF method. For comparison purpose, we consider also a Two- 
Step Method (TSM). 

The paper is organized as follows. In Section [2] we present the class of seasonal fractionally integrated 
processes with stable innovations. Section [3] addresses the ECF method and the asymptotic properties. 
Section previews the Two-Step Method. Section [5] illustrates the ECF procedure and compares it to the 
TSM in a Monte Carlo study. Section|6]concludes. 

2. Model 

2.1. Stable distributions 

In this section, we summarize the relevant facts associated with the stable distributions and refer the reader 
to Samorodnitsky and Taqqu ([23]) for a detailed statistical description. There are several ways of defining 
the stable distributions. The a-stable distribution requires four parameters for complete description: an 
index of stability a e (0, 2] also called the tail index, tail exponent or characteristic exponent, a skewness 
parameter £1 e [-1, 1], a scale parameter cr > and a location parameter /j e R. The tail exponent a 
determines the rate at which the tails of the distribution tape off. When a = 2, the Gaussian distribution 
results. When a < 2, the variance is infinite and the tails are asymptotically equivalent to Pareto law, i.e. 
they exhibit a power-law behavior. When a > 1, the mean of the distribution exists and is equal to /j. When 
the skewness parameter /3 is positive, the distribution is skewed to the right and when it is negative, it is 
skewed to the left. When (3 — 0, the distribution is symmetric about jj. As a approaches 2, loses its 
effect and the distribution approaches the Gaussian distribution regardless of The last two parameters, 
cr and are the usual scale and location parameters, i.e. cr determines the width and fi the shift of the 
mode of the density. In general, it will be convenient to define the a-stable random variables in terms of 
their characteristic functions. A random variable X is said a-stable, denoted S a fi([i, cr), if its characteristic 
function is given by 



<t>x(0 = • 



' exp [ifit - cr" \t\ a [l - ij58, tan *f 
exp{///f-cr|f|[l +;£|c;,ln|f|]} 



]} ifff^l, 



if a = 1, 



(2.1) 



where 6, 



I if t > 0, 
if t = 0, 
-1 if r < 0. 
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Using the inverse Fourier transform of the characteristic function Q>x(t), we can give a integral representa- 
tion of the probability density function as: 

i r + °° 

f(x;a,o;j3,iz) = — exp(-to)O z (f) df. (2.2) 
2tt J-oc 

Then the probability density function of a symmetric standard a-stable random variable (i.e. cr - 1, fi = 
and j6 = 0) can be expressed as: 

1 f + °° 

/ (x; a) = - exp (- \t\ a ) cos (fx) df . (2.3) 
Jo 

Formula ( 12.31 l does not have closed form expression, except in three cases (Levy, Cauchy and Gaussian 
distributions). However, it can be numerically integrated. 

In the following of this paper, we will consider the standard symmetric a-stable distribution that we denote 
SaS (i.e. the case where B = 0, fj. — and cr = 1). 

2.2. Stable ARFISMA model 

Seasonal Fractional Autoregressive Integrated Moving Average time series with symmetric c-stable (SaS) 
innovations, denoted hereafter by ARFISMA-So-S, were studied in Diongue et al. it 1 Ofl . These models ex- 
hibit long range dependence, seasonality and high variability, and are an infinite variance counterpart to the 
ARFISMA model introduced by Reisen et al. 112 ill . In this Section, we examine the definition and the basic 
characteristics of stable ARFISMA model, and we refer to Diongue et al. ifioll or Ndongo et al. lfl6ll for a 
detailed description. 

Suppose (Z,), eZ is a sequence of independently and identically distributed (i.i.d.) SaS (0 < a < 2) random 
variables with mean zero and scale parameter equals to 1 . Let B be the back shift operator and s the seasonal 
period, then the polynomials of orders p, q, P, Q are respectively defined by: 

(j> p {B) = 1 - 0iB - <p 2 B 2 <p p B p 6 g (B) = 1 + 6iB + e 2 B 2 + ■■■ + 6 q B q 

<b P {B s ) = 1 - <D. S B' 5 - <D 2 . s B 2s <&p s B Ps ®q(B s ) = 1 + ® S B" + & 2s B 2s + ■■■ + ® Qs B Qs . 

It is assumed that these polynomials have no common zeros and satisfy the conditions O(z 5 )0(z) + and 
®(z s )6(z) + for |z| = 1. Futhermore, in the above equations, (<&,)i<,<p, (4>j)\<j< P , (®k)\<k<Q and (0/)i</< ? 
are unknown parameters. 

A zero-mean process (X,), e % is said a seasonal fractionally integrated process with SaS innovations, denoted 
here by ARFISMA(p, d, q)x(P, D, Q) J -So-S, if the following equation is satisfied 

<P P (B)<I>p(B s )X t = (1 - BY d {\ - B s y D q (B)®Q(B s )Z„ (2.4) 

where the long-memory parameters d and D are fractional parameters. Notice that this model is a direct 
generalization of the ARFIMA-SaS model of Kokoszka and Taqqu [15] and it contains several particular 
cases (e.g. Diongue et al. [ 10] for more details). We assume that the following condition holds: 

\d +D\ < 1 - - and \D\ < 1 - -, with 1 < a < 2. (2.5) 

a a 
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According to Giraitis and Leipus 01211 or Reisen et al. 112 ill , one can easily show that 



(/ - B) d (I - B S ) D = Y\ [(1 - e a 'B)(\ - e- iA 'B)f 



[|] 

= - 2cosAjB + B 2 ) dj , (2.6) 

with d = ^p, d, i = D, for i = 1, . . . , [f ] - 1, = f , and Aj = ?f, for j = [§]. 

By means of the expansion 

[§] +oo 

Q(l - 2coji 7 B + B 2 ) di = Yj b M> v ^ Bi ' 

7=0 7=0 



where the coefficients v)) >Q are given by: 



bj(d,v)= J] C^vo)---^,^,^), (2.7) 

fo+-+i[|]=7 

where of = (c/o, ■ • ■ , <^[4]X v = (vo, • ■ ■ , V[»]) with Vj = cos(A/), for j = 0, . . . , [I]. The weights (C/ (J,-, v,))/ e z 
are the Gegenbauer polynomials and they can be computed using the following recursion formula: 

C (di, v,) = 1 
C\ (di, v;) = IdiVi 

Cj (d h vd = 2vi + l) C hl (d h vd - (2^ + l) C hl (dt, vd , V/ > 1. 
Hence the process defined by (12.4b can be rewritten as: 

[f] 

<S> P (B S )4> P (B)X, = ® Q (B s )0 q (B) Y\d- 2 cos(Aj)B + B 2 r d 'Z t . 

7=0 

Therefore, under conditions (12.5) Diongue et al. [10] show that the process (X t ) teZ is stationary and invert- 
ible, and the AR(oo) and MA(oo) representations are respectively given by: 

Op (B s ) (p D (B) m , . A 

= 0(1 _ 2vj B + B ) 'X, = V ex,-, (2.8) 

and 

&o (B s ) 6„(B) m . . ^ 

= -? v (1 - 2v,B + B 2 )- rf 'Z, = > c,Z f -v. (2.9) 
(Dp(BO0p(B)y h 

The coefficients (c/J >o and \cjj >Q are defined respectively by: 

+oo +oo 



and 



®q{z s )9,{z) ^ bj(d, v)z j = <b P W)<l> p b) j] c^, for \z\ < 1, 



(2.11) 



./=« 



7=0 



where the weights (nj(d, v)J >() are such that v) = bj(—d, v), with the coefficients \bj{—d, v)J >Q given 
in equation (I2.7l i. In the particular case where P = g = 00, It is easy to verify that the coefficients (c,) > 
and (cj^j >() can be computed using the following recursion formula: 



min(j,p) 



min(j,q) 



CO 



1 and Cj = 7Tj(d,v)- ^ <t>0 H (d, v) + ^ 0,-c/-,-, Vj > 1 



and 



tmn(j,q) 



CO 



1 and c ; - = bj (d, v) + V $iC H - V (<i, v) , Vj > 1 . 



The power transfer function of (X t ) teZ is given by: 

2 



*iH) = 



e q (e- iA )& Q (e- iAs ) 



<t> p (e- a )<bp (e- iAs ) 



2sing) 2sin(|) 



-2D 



, —n < A < 7T. 



(2.12) 



(2.13) 



(2.14) 



3. ECF method and asymptotic properties 

3.1. ECF method 

Because the Empirical Characteristic Function (ECF) is the Fourier transform of the empirical distribution 
function, it retains all the information in the sample but can overcome difficulties arising from the likeli- 
hood. The basic idea for the ECF method is to minimize some distance measure between the empirical 
characteristic function and the Characteristic Function (CF). It should be noted that, there are various ECF 
methods and the approach in the dependent case is not exactly the same as in the i.i.d case, because the 
dependence must be taken into account. In this Section, we summarize the approach in the dependent case 
and refer the reader to Knight and Yu 01411 or Yu IT26ll . for more details. 

Given a stationary symmetric c-stable ARFISMA process (X,) teZ defined by equation ( 12.4b . We denote by 
iff - (a, d, D, <pi,...,<p p , 6i,...,6 q , <J>i,. . .,Op, ©i,...,0q), its parameters. The preceding discussion 
in Section l2~2l motivates the choice of our parameter space ^ given by 

¥ = {<A e : cf> p + 0, 9 q +0,<b P + 0, e * 0, 6 q (z) , <p p (z) , <D P (z s ) and G (z s ) , 

have no common zeros for \z\ < 1, \d + D\ < 1 and \D\ < 1 , 1 < a < 2). 

a a 

We assume that if/o is the true value of ift and that t//o is in the interior of the compact set "S c RP+? +/>+ 2+ 3 . We 
wish to estimate tf/ from a finite realization (X[ , . . . , Xt) by the ECF method. In the dependent case, the pro- 
cedures involve moving blocks of data. Denote the moving blocks for (X u ..., X T ) as Yj = (Xj, X j+m ) , j = 
1, . . . , T - m. Thus each block has m + 1 observations. 



Expressions of these coefficients for others values of P and Q can be found in Ndongo fl7ll 
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• The CF of each block is defined by: 

c(r,iff) = E (exp \ir' YAj , where r — (fj , . . . , r m+ \)' is the transform variables . 

• The joint ECF is given by: 



1 " 

c « M = _ ^ ex P i} r '^i) ' wnere n = T 



To estimate the parameter ^ via the ECF method, one can minimize 

/„(«A) = J.-.J \c„(r)-c(r,if,)\ 2 g(r) dr, 
where g(r) being a continuous weighting function. Or equivalently, one can minimize 

I H W) = J...J |c„(r) - c(r; <A)| 2 dG(r), 
or solve the following estimation equation 

J ... J (c„(r) - c(r; if,)) W (r, i/r) dr = 0, 
where G(r) and W(r, if/) are weighting functions. 



(3.15) 



(3.16) 



(3.17) 



To use the ECF method, we have to ensure, at first, that the joint CF of the stable ARFISMA model has 
a close form. On the other hand, we need to determine the value of the size of moving blocks m and the 
weighting function to be used. Thus, we show in theorem below, that the joint CF has close form expression. 



Theorem 3.1. Let (X,), eZi be an ARFISMA-SaS process defined by equation ( 12.4l i. The joint CF of X,- m , X t -„ 
is given by: 

m+l m+l-l 



c(n,r 2 ,..., r m+ i ; iff) = exp < - ^ 



Z 

1=0 

where the coefficients (c,) are given by equation (12.1 U 



n+iCj+i 



Z Z 

1=2 h=0 



Proof. The proof is postponed to Appendix A. ■ 

Concerning the choice of the weighting function, according to Knight and Yu [ 14] or Yu ll26Tl one can obtain 
a optimal weighting function W*(r, iff) given by: 



W (r, iff) 



.dX 



J+m> 



(3.18) 



where / (Xj +m \Xj, . . ., Xj +m -ij is the conditional score function. This weight is optimal in the sense that the 
asymptotic variance of the estimator based on equation (13 . 17b can be made arbitrarily close to the Cramer- 
Rao lower bound when m is large enough. Because the conditional score of the stable ARFISMA model is 
unknown, these procedure is not feasible. Instead, in this paper, we use the sub-optimal ECF method with 
an exponential weight g{r) = exp(-r'r). The advantages of using an exponential weight are twofold. First, 
it puts more weight on the interval around the origin, consistent with the recognition that the CF contains 
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the most information around the origin and the second reason is for computational convenience. 



It is very important to recognize that when using the joint ECF, an additional choice needs to be made, 
which is that of the overlapping size of the moving blocks, m. The efficiency of the ECF estimator is very 
sensitive of the choice of m, as the moving blocks with a different m may contain different amounts of 
information in the sample. In general, there is a trade-off between large m and small m. Since the process 
in ( 12.4b does not have a Markov property, a large value of m works better than the smaller m in term of the 
asymptotic efficiency, because it retains all the information in the moving blocks of the original sequence. 
Unfortunately, it is not very clear how m affects efficiency, and we do not have an obvious guide to the 
choice of m. Ideally, m should be selected so that it minimizes the mean-square error (MSE) of the ECF 
estimator. Precisely, let ^ (m) be the ECF estimator for a given m. Then we look for m opt such that 



l opt 



1 R 2 



(3.19) 



(m) 



where R is the number of replication, ipa is the nominal value of ip and \f/j is the estimate for sample j. 
3.2. Asymptotic properties 



The asymptotic properties of the ECF estimator in the i.i.d case have been obtained by Heathcote 11311 . In 
the dependent case, Knight and Yu 1141 establish, under some regularity conditions, that the ECF estimator 
is consistent and asymptotically normally distributed. They give sufficient conditions to check the con- 
ditions listed by Newey and McFadden 11811 . namely compactness, continuity, uniform convergence, and 
identifiability, under which an extremum estimator is consistent and asymptotically normally distributed. In 
this Section, our work consist to check a set of sufficient assumptions to verify regularity conditions listed 
in Knight and Yu 1 14]. So for any fixed m, the following assumptions are necessary: 

Assumption 3.2. The parameter space *P c R/> + 4 +/>+ 2+ 3 is a compact set with i//o is in the interior of *T. 
Assumption 3.3. The coefficients \ cjj >Q in the moving average representation (12. 9t satisfy the following 



properties: 



(a) 
(b) 
(c) 



Eh 



2(a-l) 



< oo 



V 1 < a < 2, 



7=0 



E 

7=0 

OO 

E 



8 Cj 



86 



ds ]og 



and Yj 



7=0 



8 2 Cj 



8686' 



7=0 

where 6 = (d, D, <p\ 



1=0 



< oo. 



,<p p , e u ...,e q , cDi,...,(D/>, ©i,...,@ Q ). 

Assumption 3.4. J . . . J" \rj\ 6 g(r) dr < oo, V j = 1, . . . ,m + 1. 

Assumption 3.5. Let 7 (i//) = J ' . . . J \c(r; <// ) - c(r; i//)| 2 g(r) dr and I (tft) = only if ip = 



1 



Assumption 3.6. Let 3j be a cr-algebra such that \Kj, 3^} is an adapted stochastic sequence, where Kj 
K{Yf, if/) and K(x; if/) is a function defined by: 



K(x; iff) = J ■ ■ ■ J (cos(r'x) - c(r; ^ ) ) -77 < ( r. ib r) dr. 



dif/ 



We can think of 3 ; as being the cr-algebra generated by the entire current and past history of Kj. Let 
Vj = E [K \K-j, K-j-u . . .] - E \Kq\K-j-u K_ h2 , ■ . .], for 7 > 0. Assume that E [K \S-i] converges in mean 
square to as I — > 00 and Y,°jLo E [v'yVyJ < 0. 



Assumption 3.7. g (r) is a probability density function in 



Remark 3.8. 

• Assumptions l3.3l and l3.4l are technical assumptions and allow to verify the hypothesis (A2), (A3), (A5) 
and (A 6 ) of Knight and Yu (3]. 

• Assumptions 13.21 [331 13.61 and 13.71 are respectively similar to assumptions (A\), (A4), (A7) and (A§) of 
Knight and Yu B1411 . Assumption ^. 21 ensures the compactness of the parameter set. Assumption 13. 5 1 
is the identification condition whereas Assumption l3.6l provides sufficient conditions for a strong law 
of large numbers and a central limit theorem for a strictly stationary and ergodic sequence. 

The asymptotic properties of the ECF estimator based on the joint characteristic function are summarize in 
the following theorem: 



Theorem 3.9. Let tf/„ = argmin I„(if/). Suppose that Assumptions 13. 2113. 51 and Assumption 13771 hold, then 

A P- s 1 1 

if/,, — > if/Q. If, in addition Assumption l3.6l holds, then: 

Vn(^« - <Ao) N(0, B- I CAo)A(tAo)B- I ((Ao)), 

00 

where A(^ ) = Var(K{Y x ; if/ )) + 2^cov(K(Yu if/ Q ), K(Yf if/ )) 

7=2 

and B~ l (if/u) is the inverse of the matrix B(tf/o) defined by: 

#(<Ao) = J ■ ■ ■ J ^ c(r; ^0)^7 c(r; if/ )g(r) dr. 
Proof. The proof of this theorem is given in Appendix B. ■ 



4. The Two-Step Method 

In this section, the Two-Step Method (TSM) is used to estimate the parameters of ARFISMA-So-S model. 
Suppose that X\,. . . ,Xj are generated by the ARFISMA-SaS process defined by equation (12.4b . Denote 
by £ = (d, D, <pi, . . ., <p p , 9\,...,6 q , <Di , . . . , Op, ®i , . . . , &q), and if/ = (a, the parameter vector to be 
estimated. In the first step, we use the MCMC Whittle method developed in Ndongo et a /. II 1 al to estimate 
the parameter vector f. In the second phase, the MLE method given in Alvarez and Olivares U] can be 
employed in the filtered series to estimate the innovation parameter a. Therefore, the Two-Step Method can 
be described as follows. 
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4.1. The MCMC Whittle method 

This method is based on approximation of the Whittle likelihood function, using Markov Chains Monte 
Carlo (MCMC) method. Thus, the MCMC-Whittle's estimator Xw of £ is obtained by minimizing the 
following likelihood function: 

1 N 1 

L W (X, £) = -Y___ (4.20) 
N j-J h x (Aj,0 

where N is taken large enough from the strong law of large number and hx(A, is the power transfer 
function of the process (X r ) teZ generating the data, and is defined by equation (12.141 1. In equation (14.201 ). the 
sequence Ai,...,An is generated using a Metropolis-Hastings algorithm (see Ndongo et al.\\6] for more 
details). 



4.2. The EML method 

In this second step, we firstly use the estimation obtained by MCMC- Whittle method to calculate the filtered 
series (Z,) te z, thanks to the autoregressive representation (12.8) . However, given the observation X\, . . . ,Xj 
the innovation (Z,) l=l T cannot be directly computed, since an infinite sample would be needed. Neverthe- 
less, they may be estimated by: 

T-l 

Zt{tw) = J]Zj(iw)Xt-j, t=l,...,T, 
j=o 

where the coefficients (cj(£w)) >0 arc defined by equation ( 12.10b . Secondly, we estimate the innovation 
parameter a given the observations Z\, . . . , Zj by using a MLE procedure developed in Alvarez and Olivares 
yj]. Therefore, a is obtained by maximizing the following log-likelihood function: 

T 

1(a) = J] log (f(Z t , a)), 
where f(Z,,a) is the density probability function for the random variable (Z,),<=z- 



5. Monte Carlo experiments 

The Monte Carlo study is designed to check variability of the ECF procedure in comparison with the Two- 
Step Method. The simulation results give the average values, the root mean square error (RMSE) and the 
mean absolute error (MAE) of the estimation procedures based on 1500 replications. The number of obser- 
vations is set at T = 1500 and the seasonal period is fixed at s = 4. All calculations were carried out using 
an R programming environment (see N20ll ) on a Intel(R) Pentium Dual-Core Processor T2050, 1.60 GHz (2 
CPUs) computer. 

We generate different ARFISMA sequences with SaS innovations. We adopt the method developed by 
Stoev and Taqqu B2411 for FARIMA time series with SaS innovations. We approximate the path X,, t = 
1 , . . . , T by the truncation moving average 

M 

Xt = YjCjZ H , (5.21) 

7=0 
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where the non random constants (cj)j>o are defined by equation (12.1 11 , and M is the truncation parameter 
and fixed in this study to 5000. Indeed, in order to represent well the long-term dependence behaviour, 
Stoev and Taqqu J24ll suggest to use large values for the truncation parameter. They show also that the 
choice of M greater or comparable to the sample T work well in most case. The SaS innovations are 
obtained through a version of the Chambers et al. |5J algorithm described in Section 1.7 of Samorodnitsky 
and Taqqu lEjI . All of the simulations involve one of the following models summarized in Table Q] 



Parameters a d D <p 6 

Model 1 L6 Oil O20 

Model 2 1.6 0.15 0.20 0.60 

Model 3 1.6 0.15 0.20 0.40 

Model 4 1.6 0.15 0.20 0.60 0.40 



Table 1: Data generating processes. 



For each model, we choose, at first, several values for the size of the moving blocks m (m — 1 to m = 8) to 
examine the effect of m on the estimates. Then, we select the optimal m for each model using the criterion 
described by equation ( 13.191 1, and we compare the results with the Two-Step Method. On the other hand, 
we vary the values of short memory parameters, to study their effects on the estimation of long memory 
parameters. We show here how to compute the ECF estimator in practice for m = 1, and the procedure is the 
same for other values of m. Note that with m = 1, the characteristic function and the empirical characteristic 
function are respectively given by: 



c(ri,r 2 ;t//) = exp 



-Z 

7=0 



ncj + r 2 Cj+i\ -\r 2 \' 



and 



1 " 

c„(ri,r 2 ) = ~ Zj ex P( l ' r i- x 'j + ir 2 x j+i) . n = T-l. 



(5.22) 



(5.23) 



In this paper, we use the ECF method with an exponential weighting function, then the procedure is to 
choose iff to minimize 



JJ' 



\c„(n, r 2 ) - c(r u r 2 ; i/0| 2 exp(-r 2 - r 2 2 ) dr x dr 2 , (5.24) 

where c(ri,r2',i//) and c n (r\,r2) are given by ( 15.221 1 and (15.231 1. The implementation of the ECF method 
essentially requires minimizing (15.24b . which involves double integrals. Consequently, we use a Monte 
Carlo Integration method to evaluate the multiple integral ( 15.24b . 



5.7. Empirical choice of the size of the moving blocks 

In this Section, we show the effects of the size of moving blocks m and, then we determine the optimal 
value of m suitable for each model presented in Table Q] The simulation results are given in Tables [2] [3] |4] 
and|5] 

• In the first experiment, a symmetric <r-stable ARFISMA(0, d, 0)(0, D, 0) model is considered (i.e the 
model without short memory parameters), and the results are summarized in Table [2] We see that, 
for all the values of m, the ECF method performs very well as the RMSE and the MAE are small. 
However, this performance varies with the value of m. Indeed, for the long-memory parameters d 
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and D, the best estimates are obtained for m = 1 whereas for the innovation parameter a the value 
m — 4 is the best one. Nevertheless, in terms of RMSE and MAE the difference between m = 1 and 
m — 4 for the parameter a is relatively small (respectively 1.11% and 1.18%). Therefore, we consider 
the value m — 1 as the optimal value of the size of moving blocks to estimate simultaneously all the 
parameters of Model 1 . 

• In the second experiment, we introduce in the previous model a short memory part, and the results are 
displayed in Tables |3]|4] and [5] Firstly, in these tables we observe that the RMSE and the MAE of the 
estimates are larger for all values of m compared to those obtained in Table [2] Consequently, the AR 
and MA components can also be a source of the bias of the estimators. This phenomenon is already 
observed in the literature (e.g. Boutahar et al. [2], Diongue and Guegan |@]). On the other hand, the 
results are very bad for m = 1, 2, 3 and 4, and we observe a significant improvement from m = 5. 
In general, we remark that the value m - 6 seems to give better performance. Hence, we consider 
the value m — 6 as the optimal value of the size of moving blocks to estimate simultaneously all the 
parameters of the models with short memory components. 

In view of the simulation results, we can conclude that the efficiency of the ECF method depend on the 
choice of the size of moving blocks. Hence, we need to determine the optimal value of m. Thus, for the 
stable ARFISMA(0, d, 0)(0, D, 0) model, i.e. without short memory components, the value m = 1 is the 
optimal value of m whereas when the model has a short memory part the value m — 6 is the best one. In the 
Section below, the estimates with these optimal values of m are compared to those obtained by the Two-Step 
Method (TSM). 





m — 1 






m — 2 






Statistics 


a 


d 


D 


~a 


d 


D 


Mean 


1.615 


0.132 


0.170 


1.594 


0.130 


0.155 


RMSE 


0.015 


0.017 


0.029 


0.006 


0.019 


0.044 


MAE 


0.015 


0.017 


0.029 


0.005 


0.019 


0.044 




m — 3 






m — 4 








ce 


d 


D 


7r 


d 


D 


Mean 


1.573 


0.128 


0.136 


1.596 


0.128 


0.145 


RMSE 


0.026 


0.021 


0.063 


0.004 


0.021 


0.054 


MAE 


0.026 


0.021 


0.063 


0.004 


0.021 


0.054 




m — 5 






m — 6 








~a 


d 


D 


a 


d 


D 


Mean 


1.591 


0.131 


0.144 


1.593 


0.129 


0.144 


RMSE 


0.009 


0.018 


0.055 


0.007 


0.020 


0.055 


MAE 


0.008 


0.018 


0.055 


0.006 


0.020 


0.055 




m = 7 






m = 8 








c? 


d 


D 


a 


d 


D 


Mean 


1.596 


0.127 


0.145 


1.601 


0.125 


0.148 


RMSE 


0.005 


0.022 


0.054 


0.006 


0.024 


0.051 


MAE 


0.004 


0.022 


0.054 


0.004 


0.024 


0.051 



Table 2: ECF estimation for model 1, for m = 1 to m = 8. 
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m — 1 








m = 2 








Statistics 


a 


d 


D 


Y 


~a 


d 


D 


d> 

T 


Mean 


1.726 


0.217 


0.233 


0.464 


1.683 


0.125 


0.252 


0.633 


RMSE 


0.288 


0.208 


0.263 


0.366 


0.171 


0.117 


0.132 


0.192 


MAE 


0.168 


0.174 


0.147 


0.292 


0.122 


0.100 


0.089 


0.162 




m — 3 








m — 4 










a 


d 


D 




a 


d 


D 


d> 

T 


Mean 


1.663 


0.125 


0.243 


0.641 


1.650 


0.128 


0.235 


0.644 


RMSE 


0.135 


0.109 


0.101 


0.167 


0.120 


0.106 


0.087 


0.154 


MAE 


0.100 


0.095 


0.073 


0.140 


0.085 


0.093 


0.063 


0.130 




m — 5 








m = 6 










<? 


d 


D 


4> 


a 


d 


D 


<P 


Mean 


1.634 


0.160 


0.228 


0.606 


1.624 


0.169 


0.222 


0.591 


RMSE 


0.098 


0.097 


0.070 


0.119 


0.082 


0.099 


0.057 


0.103 


MAE 


0.068 


0.078 


0.047 


0.089 


0.061 


0.079 


0.040 


0.083 




m = 7 








m = 8 










~a 


d 


D 






d 


D 


<P 


Mean 


1.571 


0.250 


0.187 


0.394 


1.568 


0.248 


0.189 


0.395 


RMSE 


0.092 


0.135 


0.059 


0.258 


0.093 


0.132 


0.056 


0.247 


MAE 


0.065 


0.114 


0.043 


0.219 


0.067 


0.110 


0.041 


0.210 



Table 3: ECF estimation for model 2, for m = 1 to m = 8. 





m = 1 








m — 2 








Statistics 


a 


d 


D 


e 


a 


d 


D 


8 


Mean 


1.694 


0.177 


0.250 


0.339 


1.672 


0.165 


0.240 


0.399 


RMSE 


0.209 


0.064 


0.130 


0.153 


0.156 


0.056 


0.094 


0.059 


MAE 


0.133 


0.046 


0.077 


0.090 


0.112 


0.041 


0.067 


0.041 




m — 3 








m = 4 










a 


d 


D 




<? 


d 


D 


8 


Mean 


1.659 


0.169 


0.239 


0.394 


1.637 


0.170 


0.230 


0.398 


RMSE 


0.135 


0.057 


0.084 


0.052 


0.102 


0.059 


0.070 


0.054 


MAE 


0.097 


0.042 


0.059 


0.039 


0.074 


0.044 


0.047 


0.041 




m = 5 








m = 6 










a 


d 


D 


8 




d 


D 


8 


Mean 


1.630 


0.170 


0.225 


0.398 


1.534 


0.126 


0.143 


0.343 


RMSE 


0.095 


0.064 


0.066 


0.059 


0.066 


0.023 


0.056 


0.057 


MAE 


0.067 


0.046 


0.043 


0.044 


0.065 


0.023 


0.056 


0.056 




m = 1 








m = 8 










a 


d 


D 


6 




d 


D 


8 


Mean 


1.622 


0.171 


0.223 


0.396 


1.619 


0.170 


0.222 


0.396 


RMSE 


0.090 


0.071 


0.070 


0.067 


0.083 


0.069 


0.075 


0.068 


MAE 


0.063 


0.049 


0.043 


0.051 


0.061 


0.050 


0.044 


0.053 



Table 4: ECF estimation for model 3, for m = 1 to m = 8. 
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m — 1 m = 2 



Statistics 


(Y 


d 


D 




6 


CY 


d 


D 


d> 


e 


Mean 


1.488 


0.356 


0.0008 


0.002 


0.331 


1.508 


0.275 


0.026 


0.276 


0.411 


RMSE 


0.111 


0.206 


0.199 


0.597 


0.069 


0.091 


0.125 


0.173 


0.323 


0.012 


MAE 


0.111 


0.206 


0.199 


0.597 


0.068 


0.091 


0.125 


0.173 


0.323 


0.011 




m — 3 










m — 4 










Statistics 


(T 


d 


D 


4> 


e 




d 


D 


4> 


e 


Mean 


1.513 


0.262 


0.045 


0.288 


0.402 


1.514 


0.222 


0.089 


0.349 


0.365 


RMSE 


0.086 


0.112 


0.155 


0.312 


0.005 


0.085 


0.073 


0.110 


0.250 


0.034 


MAE 


0.086 


0.112 


0.154 


0.311 


0.004 


0.085 


0.072 


0.110 


0.250 


0.034 




m — 5 










m — 6 










Statistics 


~a 


d 


D 


<P 


9 


~a 


d 


D 




6 


Mean 


1.520 


0.178 


0.118 


0.451 


0.351 


1.513 


0.147 


0.133 


0.510 


0.334 


RMSE 


0.080 


0.031 


0.081 


0.151 


0.049 


0.087 


0.017 


0.066 


0.095 


0.066 


MAE 


0.079 


0.028 


0.081 


0.148 


0.048 


0.086 


0.013 


0.066 


0.089 


0.065 




m = l 










m = 8 










Statistics 




d 


D 


<P 





<? 


d 


D 




6 


Mean 


1.492 


0.135 


0.139 


0.525 


0.305 


1.456 


0.132 


0.142 


0.503 


0.268 


RMSE 


0.109 


0.028 


0.062 


0.089 


0.096 


0.148 


0.041 


0.060 


0.124 


0.135 


MAE 


0.107 


0.022 


0.060 


0.077 


0.094 


0.143 


0.032 


0.057 


0.105 


0.131 



Table 5: ECF estimation for model 4, for m = 1 to m = 8. 



5.2. Monte Carlo study comparing ECF procedure and Two-Step Method 

In this section, we compare the finite sample performance of the ECF method to the Two-Step Method 
discussed in Section|4] The simulation results are summarized in Tables |6]|7] [8] and [9] 

• When we consider Table|6]displaying the estimation results for ARFISMA(0, d, 0)X(0, D, 0)j-SaS 
model, we see that all methods perform well, as the RMSE and MAE are in most cases small. The 
estimation parameters from the TSM are better than those obtained by the ECF approach. 

• Tables [7] [8] and [9] display the estimation results, when there are long memory and short memory 
components simultaneously. From these tables, we observe that the estimation results are acceptable 
for all methods, in most cases. However, it can be remarked that the estimation of the parameters is 
disturbed. Indeed, the RMSE and MAE obtained from these tables are larger than those presented 
in Table [6] Hence, this shows that the impact of the short memory parameters in the estimates. 
Comparing the two methods, it can be seen that better estimates are obtained, in general, from the 
ECF procedure than the TSM. Moreover, an advantage of the ECF method is that, as a consistent 
procedure it estimates all of the parameters simultaneously. 



Statistics 


ECF (m= 


1) 






TSM 




a 


d 


D 


a 


d 


D 


Mean 


1.6158 


0.1321 


0.1707 


1.6105 


0.1565 


0.1844 


RMSE 


0.0158 


0.0178 


0.0292 


0.0108 


0.0132 


0.0193 


MAE 


0.0158 


0.0178 


0.0292 


0.0105 


0.0106 


0.0164 



Table 6: Monte Carlo study to compare ECF procedure with TSM for model 1. 
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ECF (m=6) 






TSM 






Statistics 


a d 


D 4> 


a 


d 


D 


<P 


Mean 

RMSE 

MAE 


1.624 0.169 
0.082 0.099 
0.061 0.079 


0.222 0.591 
0.057 0.103 
0.040 0.083 


1.603 
0.008 
0.006 


0.164 
0.164 
0.077 


0.200 
0.079 
0.026 


0.587 
0.104 
0.073 


Table 7: Monte Carlo study to compare ECF procedure with TSM for model 2. 




ECF (m=6) 






TSM 






Statistics 


a d 


d e 


a 


d 


D 


8 


Mean 

RMSE 

MAE 


1.534 0.126 
0.066 0.023 
0.065 0.023 


0.143 0.343 
0.056 0.057 
0.056 0.056 


1.256 
0.344 
0.343 


0.157 
0.064 
0.026 


0.204 
0.068 
0.022 


0.395 
0.048 
0.029 



Table 8: Monte Carlo study to compare ECF procedure with TSM for model 3. 







ECF (m 


=6) 










TSM 






Statistics 


a 


d 


D 


<P 





a 


d 


D 




e 


Mean 

RMSE 

MAE 


1.513 
0.087 
0.086 


0.147 
0.017 
0.013 


0.133 
0.066 
0.066 


0.510 
0.095 
0.089 


0.334 
0.066 
0.065 


1.286 
0.316 
0.313 


0.146 
0.050 
0.040 


0.182 
0.025 
0.021 


0.610 
0.062 
0.049 


0.172 
0.233 
0.227 



Table 9: Monte Carlo study to compare ECF with TSM for model 4. 



6. Conclusion 

In this article, we studied the ECF method developed by Knight and Yu [14]. The method is illustrated 
to estimate simultaneously the parameters of stable ARFISMA models introduced by Diongue et al. [10]. 
Under some conditions, the resulting estimators are shown to be consistent and asymptotically normally 
distributed. For comparison purpose, we have also consider a Two-Step Method (TSM) composed by the 
MCMC Whittle procedure developed in Ndongo et al. 11611 and the MLE approach introduced in Alvarez 
and Olivares HI. Finite sample behaviours of these methods were studied through Monte Carlo simula- 
tions. It is found in general, that the ECF method is better than the TSM, particularly, when short memory 
components are present in the model. The simulation results show also the impact of the size of moving 
blocks and the short memory parameters in the estimate. 

It is very important to remark that there are two difficulties with the implementation of the characteristic 
function-based estimators (the use of optimal weighting function and large set of moment conditions) as 
noted in Carrasco et al. ([6]). They solved the two problems in the framework of the GMM with continuum 
of moment conditions. However, the models considered here is non-Markovian, and then the conditional 
Characteristic Function (CCF) is unknown and difficult to estimate. Hence, the GMM based on the CCF 
is not feasible. On the other hand, the Joint Characteristic Function (JCF) of stable ARFISMA models is 
known but the GMM -JCF will not be efficient, since the process is not Markovian (see Carrasco et al. (|Ql) 
for more details). Thus, it will be interesting to compare the sub-optimal ECF method with an exponential 



14 



weighting function and the GMM based on the JCF, through Monto Carlo simulations. This problem will 
be examined in a forthcoming paper. 
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Appendix A: proof of Theorem [37T1 



c (ri, . . . , r m+ i ; i/0 = E [exp {iriX,- m + ir 2 X, 
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Using equation (12.9b . we can write 

c(r u . ..,r m+1 ;t//) = E 
Developing and regrouping the terms in Z t - m -j for j = to oo, we obtain 



exp I ir\ ^ c j z t-m-j + in ^ 

I 7=0 7=0 7=0 



c (r; ifr) = E 



E 



exp 



' y//! r M c i+i) z t-m-j + i 2^ r h+1 c h Z,- m+ i +■■■ + ir m+l c Z, 



7=0 /=0 



//=() 
m+1 m+l-l 



exp j i r M Cj+i)Zf- m -j + i ^\ ^ r M c h )Z h 

I j=o /=o 



■m-l+/ | 



7=2 /i=0 



Using the independence of the Z,, we have 
c (r; i/>) = E 



exp j z J^Zj r M Cj+i)Zt- m -j 
7=0 /=o 



m+1 m+1-/ 



exp j z 2]( 2^ ni+iCh)Z,- 



m- 1 +i 



/=2 A=0 



Let 



A = E 



/>• - E 



exp | i ^(Jj r M Cj+i)Z t - m -j 
{ 7=0 /=o 

f m+1 m+l-l 

exp ^ r h+l c h )Z t . 

{ 1=2 h=0 



m-l+l { 



Calculate A and B 
B = E 



= E 



I m+1 m+l-l \ 
*2^Zj r h+iCh)Z,- m -\+i\ 
1=2 h=0 j 
m+1 f m+l-l 

j~[exp j;'( ^ r h+ ic h )Z t 

1=2 { h=0 

I ( m+1-/ 

| | - exp | rh+iCi,)Z, 



't-m-l+l 



-m-l+l } 



for indenpence of the Z, 



m+1 (m+l-l \ 

rf,+iCh 

1=2 V h=0 



m+l-l 



- n*- z 

1=2 \ h=0 
m+ 1 ( r 

n-p - z 

V A=0 
f m+1 m+l~ 

-z z 



/=2 



= exp< 



rh+iCh 



1=2 h=0 



For calculation of A, let F = r i+i c j+i) z t-m-j and remark that A = <I>y(l), where ®y is the CF of Y. 

>=o /=o 



The sequence {Z,} is i.i.d symmetric a-stable and ^ 2 r '+ lC 7+' 

7=0 '=o 
16 



oo, then according to Cline (01) or 



Cline and Brockwell ([8]), we have Y is also symmetric <?-stable and 



Z Yj r i +l c j+ i 

V./=o /=o 



l/ar 



Zi. 



As consequently 



<D f (m) = E 



= E 



7 



exp j z'k Z*Z r l+iCj+i)Z,- m -j 



j=o 1=0 



exp < ;m 



Z Xi r ' +ic > +/ 

V;=0 (=0 



1/a 



Z Z r ' +ic ^ + ' 

V;=0 (=0 



l/a\ 



= expj-|«rZ Z r/+lC 

7=0 /=0 



-j+l 



Hence A = <D y (1) = exp I - ^ ^ r, +I c i+ , 

1 7=0 /=0 

Finally 

c(r;ip) = Ax5 

J 00 m 

= exp rZ Z r, + ic >+' f xe H~Z Z 



m+l 



/=2 



m+l-/ 



/i=() 



= ex pj-Z Z r,+ic > +/ 

I ;=0 1=0 



m+l-/ 



I I 

/=2 fc=0 



Oh-/ca 



Appendix B: proof of Theorem l3.9l 

In this Section, we present the proof of Theorem 13 . 9 1 which consists, under Assumptions l3.3l and l3.4l to 
verify the hypothesis of Theorem 2. 1 of Knight and Yu B1411 . namely Assumptions (A2), (A3), (A5) and (A(,). 

(A2): With probability one, I„{ifi) is twice continuously differentiable under the integral sign with respect 
to ip over *P. 



Let /(r; if/) - \c n (r) - c(r; ip)\ g(r), then we have: 



Olp 



2g{r)\c n {r)-c{r,iff)\ 



d c(r; ip) 



< 2g(r)(\c„(r)\ + \c(r,if/)\) 
d c(r; if/) 



dip 
d c(r; if/) 



dif/ 



< 4g(r) 



dip 
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Defining 



U{r;xp) = 2_jYj n+lC j +l 
j=0 1=0 



m+1 m+l—l 

+ ^ 2^1 n,+lCh 

1=2 h=0 



we can rewrite the previous inequality as: 
8 



dip 



fin iff) 



< 4g(r)exp(~E/(r;#)) 
8 U(r; iff) 



8 U(r; i/O 



dip 



< 4g(r) 



dip 



• Let us show that I„(ip) is of class C 1 with respect to 5 = (d, D, 

(01,... ,^),fl=(fl lf ... ,«,),* = (*!,..., Op) and = (©!,..., e ). 



Let5(r;0) = ^ K,(r; 0), with Uj(r, ip) = ^r M c j+ i 
7=0 /=o 



Then we get: 



3 Uj(r; ip) 



88 



/ in 



n+iCj+i 



8 cj+i 



< 2 



1=0 



<9<5 



^In+iCj+il ^k/+iflj+/l , withflj = — -. 

u=o / w=o / 

Applying Holder's inequality to the previous relationship, we obtain: 



3 Uj{r; ip) 



86 



< 2 



< 2 



m / m / hi \ 



V/=0 
/ m 



V /=() 



V/=0 



2 I J]\c j+ ,\ a J]\a j+l \° 



V/=0 



/ V 1=0 



V/=0 



■ There exists Iq such that \Cj+i\ a < \cj + i \ a , for all I = 0, . . . ,m. 



V/=0 



^(m+l) 1 "!^!"- 1 . 



There exists such that |a ;+ /| u ' < \cij+i { I", for all / = 0, . . . ,m. 



8, <t>, ©), where 



(6.25) 



v /=() 



< (m + l)«|a 7 - + ;j|. 



Using expression found in ( 16.25b and ( I6.261 l. we get 

< 2(m + 1) 



8 uj(r, ip) 



86 



< (m+1) 



v;=o ^ 

1^ 



V/=0 



Wy, With Wj = \c j+k \ 2( > a -V 



(6.26) 



18 



According to Assumption 13. 31 the series 2 Wj is convergent. Consequently, we get U(r; if/) is differ- 
entiable with respect to 5 and we can write: 

a-1 



d U(r; \fr) 



35 



88 

< (m+1) 



m+ 1 m+l-l 

3S(r\ip) V V 

1=2 h=0 
( m \a co 

Z |r,+i1 Z w > 

1=0 > j=0 

+1-1 \°- ] 
^ k/ J+ ;Ql 2j k/i+za/il 



m+l-i 



/!=() 



/=2 



/;)!+!-/ 



V /!=0 



< M (m + 1) 



V /!=0 
\« m+1 



Z |r,+i1 +a Z 

/=0 / 1=2 



(m+l-l 



\ a ~ x (m+l-l 

^ k/ !+ /Qil ^ kft+z«ftl 



V /i=o 



V /i=0 



m+1 



Let A = g 2^ 

/=2 



/m+l-i 



V /i=o 



(m+\-l 



2^ k/ !+ /Qil 2j k/i+ifl/il 



V /i=0 



Applying Holder's inequality to A, we have: 



/m+1-/ V*- 1 
, /i=0 

(m+l-l V'"' 1 
2 k/,+/| 



V /i=0 



A < 2 J] 

/=2 
m+1 

/=2 
m+1 

- 2 Z Z |r ' !+ '' Z w 

1=2 IA /i=0 

m+1 f m+l-l 

< 2M'£ J> 



/=2 ^ 

m+1 r (m+l-l \ a t +oo ^1/2 

«/il 

/=2 IA h=0 ) U=o 

m+1 (m+l-l \ a 



'm+l-l 
, h=0 


1/2 


'm+1-/ 

, /i=0 , 


'm+l-l , 
, /i=0 


'm+1-/ > 

2 l«// 
. A=0 


1/2- 



1/2 



V h=0 



using Assumption l3.3 



So we get: 

5 U(r; iff) 



35 



m+1 (m+l-l \ a 

Yj\n+i\ +2M'J] J] \r m \ 



\i=o 



< M(m+1) 

< M (m + l)Hi (r) +2M'Y//, (r; /) 



' /=2 V ft=0 

m+1 



/=2 



where 



Z"« 



\i=o 



Z ' r/+i 



\2 



Xkwl > i 



/=o 



]T|r, +1 |, si |]k,+i| < l. 



/=o 



/=o 



(6.27) 
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and 



( (m+\-l Y 

Z low I 



/m+l-i \ a 



2 inwi 



V /i=0 



<H l (r,l) = 



V A=0 

m+1-/ 



m+1-/ 

si ^ |r A+ /| > 1 

h=0 
m+1-/ 



2 l r >wl, Si X ' r/,+ '' " L 

h=0 h=0 



m+1 



Letting Ti(r) = M (m + l)#i(r) + 2 M' V //i(r; 0, we get: 



/=2 



<95 



<4Ti(r)g(r). 



Using Assumption l3.4l it is easy to show that: 

J... J r 1 (r)g(r) fl fr< £ x,. 

Hence I„(tf/) is of class C 1 with respect to 5. 

• Let us show that I„(tfr) is of class C 1 with respect to a. 

Using the same approach as before, we show that U(r; if/) is differentiable with respect to a. Thus 
we get: 



d U{r; iff) 



da 



o o m+1 m+1-/ 



/=2 ft=0 
m+1 m+1-/ 



= T -5 ^ ^ + Z\j Z r '' + ' Ch ln Z r/,+/C/, 

a 1=2 h=0 h=0 

Applying the Holder's inequality and using the approach as before, we have: 



m+1-/ 



8 U(r; iff) 



da 



+oo m+1 
2 



< (m+l) 2 // 2 (r)^|c y+/o | 2 + 2 



7=0 1=2 
m+1 



m+1-/ 

Z 

/i=0 



< M(m+l) 2 // 2 (r) + ^// 2 (r;0:=r 2 (r) 



/=2 



where 



H 2 (r) = 



2>> 

1=1) 



\3 



^ m \ 2 
V/=0 



£k/+il > i 

1=0 
111 

£k/ + il < l, 
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and 



( fm+l-l \ 



H 2 (r,l) = 



Thus we obtain: 



z 

k A=0 

•«+ 1- 

z 



KA+/I 



m+l-( 



KA+/I 



^ kA+/l > i 

A=0 
i+l-Z 

X Inwl < i. 



A=0 



5/(r;iZV) 



< 4 r 2 (r) g(r), 

and then using Assumption 13. 41 we have . . ^F2(r) g(r) t/r < oo, and hence 7„(iA) is of class C 1 
with respect to a. 

Now we will show that I n (ip) is of class C 2 with respect to ip. Thus we get: 



jf(n<P) 



2g(r) 



3 c(r; <A) 



d c(r; (fr) 



<9^ 

2^r)|c„(r)-c(r,^)| 
d c(r; ^) 



dip' 

d 2 c(r; ip) 



< 2g(r) 

< 2g(r) 
+ 4 g(r) 

< 6g(r) 



5 C/(r; <A) 



dipdip' 
d c(r; ip) 



34) 
d U{r; iff) 



dip' 
d U(r, ip) 



+ 4g(r) 



a 2 c(r; i/r) 



dip 
d U(r, ip) 



dip' 
d U{r; ip) 



dip 



dip' 
d U(r; iff) 



d 2 U(r; ip) 



dip' 



+ 4 g(r) 



dipdip' 

d 2 U(r; ip) 



dipdip' 



• Let us show that I„(ip) is of class C 2 with respect to a. 

According to the previous calculations, we have: 



d 2 



dada 



-,f{r,ip) 



<6g(r)[r 2 (r)] z +4g(r) 



d 2 U(r, iff) 



dada' 



Using the same remarks before, we can show that U (r; ip) is twice difFerentiable with respect to a, 
and we have: 



d 2 U{r; ip) 




dada' 





d 2 „, , d ^ ^ 



dada' 



1=2 



m+l-l 



rh+iCh 



A=0 



m+1 



< M (m + l) 3 ff 3 (r) + ^ 



/=2 
m+1 



z 

A=0 



r l+h c h 



a' +2 



< M(m+ l) 3 ff 3 (r) + Ytf 3 (r;Z) := r 3 (r), 



/=2 
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where 



H 3 (r) = 



v /=() 



\3 



V /=() 



£k/ + il > i 

/=0 

2k/ + il < l. 



and 



H 3 (r,l) 



(m+\-l \ 



/ J= 

11+1-/ 
E l r ' 1+ 'l 



m+1-/ 



^ A=0 



si 2^ l r ft+/l > 1 
h=0 

m+l-Z 

si 2^ ^ L 



//=() 



Thus 



According to Assumption l3.4l we get: 



<[6r 2 2 (r) + 4T 3 (r)]g(r). 



J... J [6 T\{r) + 4 r 3 (r)] g(r) dr < oo 



and then I n (ip) is of class C with respect to a. 

• Let us show that I n (tfr) is of class C 2 with respect to 5. 

Using the same approach as before we have: 



<95d5 



< 6 #(/-)!>) + 4 #(r) 



8 2 U(r; iff) 



8686' 



Since U (r; iff) is twice differentiable with respect to 6, we get: 



8 2 U{r; iff) 



8686' 



8 S(r,tff) a v-i y^ 1 



3 2 m j(r; <//) 



./=(> 

m+1 



1=2 

Hl+l 



//=() 



m+1-/ 

E r " 



//=o 



35 



Z=2 



m+1-/ 

i 

A=0 



Z=2 



m+1-/ 



E r,+ " c " e 



ft=0 
a-2 



a-i 



m+1-/ 



O+A 



3 2 c A 



/)=() 



8686' 



m+1— I ^ m+1-/ ^ 



h=0 



h=0 



M (m + 1) 



2 

V/=0 



r/+i 



+ #i +fi 2 
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Using the same approach as before, we can compute B\ and Bi as follow: 

a-l 



m+1 



B 



- «e 

1=2 



m+1-/ 



ft=0 

m+1 /m+1-/ "\ 
i+i 



m+1-/ 



/i=0 



< J#, £ £ |r ft . 



/=2 V /i=0 



and 



m+1 /"m+1-/ "\ 
"/,+/ 



B 2 <M 1 Y j \r h . 



1=2 V h=0 



Finally we get: 

8 2 U{r; iff) 



8585' 



M (m + 1) 



M (m + 1) 



I m \° 



m+1 I m+1-/ \ 



^] r, +I + (Mi + Mi) Y ^ k/ I+ /l 



V/=0 



Z=2 V h=0 



m+1 



tfi(r) + (Mi + Mi) ^ //Kr; Q = r,(r), 



/=2 



and then 



8585 



-,f(r^) 



<6g(r)r\{ r ) + Ag( r )r x { r ). 



According to Assumption l3.4l we have J" . . . J" (6r 2 (r) + 4Fj (r)) g(r)dr < oo. Thus I„(tff) is of class 
C 2 with respect to 5. 

In short, I n (i//) is twice continuously differentiable under the integral sign with respect to i//. 
(A3): The sequence (X,) leZ is strictly stationary and ergodic. 

According to Diongue et al. II 1 Ofl - the process (X ( ) (eZ is strictly stationary and has a unique moving 
average representation given by: 

00 

= / \ c j Z ,-j, 
7=0 

where the sequence (Z,) feZ is independently and identically distributed and E/Lo | c j|° < °°> f° r a ^ 
1 < a < 2. Hence (X,) teZ is an ergodic process. Consequently, the sequence (X,) teZ is strictly 
stationary and ergodic. 
(A5): K(x; i[/) is a measurable function of x for all i/r and bounded, where 



K(x; i/r) 



J". . . J"(cos(r'x) - c(r; c(r: i;>)t,'(/\) </r. 



We get: 



(cos(r^) - c(r; </0) — c(r; <%(r) 
01^ 



< 2g(r) 

< 2g(r) 



dtfi 
8_ 
dtfi 



c(r; iff) 



Uir-if,) 



23 



Using the previous remarks, we obtain: 



<T(r) 



r,(r), si 0- = (d, D, <f>, 0, <D, 0) 
T 2 (r), si if/ = a. 



Thus 



(cos(r'x) - c(r; if/)) — c(r; ^)g(r) 



< 2r(r)g(r), V^e¥. 



Using Assumption 13. 41 we have ^ • • • J" F(r)g(r) c/r < oo, and hence K(x; if/) is a measurable func- 
tion of x for all ^ and bounded. 

r r 8 8 d 2 

(A 6 ): B(if/ ) = J ... J — c(r;if/ ) — c(r; if/ )g(r) dr is nonsingular and ^ g c(r;if/) is uniformly 
bounded by a g-integrable function over *P. 



d d 
TT c ( r '^o)T-7 c(r,if/ ) 
oif/ oif/' 



We get: 



According to Assumption l3.4l we have 

a 



— U(r;if/ )— U(r;if/ ) 



< r 2 (r) 



/■■■/ 



r (r)g(r) dr < oo, and then 



.. , c ( r >>Po)-r-- c{r; if/ Q ) 
Oif/ oif/' 



g(r) dr < oo. 



Hence fi(t^o) is nonsingular. On the other hand, remark that: 



d 2 c(r; if/) 



dif/dtfY 



Now we have: 

d 2 U(r; if/) 



d U(r; if/) 



dijj 



< r\r) + 



d U(r; if/) 



dif/' 
d 2 U(r; if/) 



d 2 U(r; if/) 



dif/dif/' 



dif/dif/' 



dif/dif/' 
This implies that: 



< T'(r) 



Fi(r), siif/= (d, D, (f>, 9, <D, 0) 
r^ir), si if/ = a. 



d 2 c(r; if/) 



dif/dif/' 



< T 2 (r) + r», 



and using Assumption l3.4l we get: 



Hence 



d 2 c(r; if/) 



dif/dif/' 



J...J (T 2 (r)+T'(r))g(r)dr<oo. 



is uniformly bounded by a g-integrable function over *P. 



Thus all conditions are verified, then we can apply the Theorem 2.1 of Knight and Yu 111411 to achieve the 
proof. ■ 

24 



